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ABSTRACT 

We use a model-independent technique to investigate the nature of the 2-15 keV X-ray spec- 
tral variability in four Seyfert galaxies and distinguish between spectral pivoting and the two- 
component model for spectral variability. Our analysis reveals conclusively that the softening 
of the X-ray continuum with increasing flux in MCG— 6-30-15 and NGC 3516 is a result of 
summing two spectral components: a soft varying component (SVC) with spectral shape in- 
dependent of flux and a constant hard component (HCC). In contrast, the spectral variability 
in NGC 405 1 can be well described by simple pivoting of one component, together with an 
additional hard constant component. The spectral variability model for NGC 5506 is ambigu- 
ous, due to the smaller range of fluxes sampled by the data. We investigate the shape of the 
hard spectral component in MCG— 6-30-15 and find that it appears similar to a pure reflection 
spectrum, but requires a large reflected fraction (R > 3). We briefly discuss physical interpre- 
tations of the different modes of spectral variability. 
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1 INTRODUCTION 

Rapid X-ray variability is common in Seyfert galaxies (Green et al., 
1993) and suggests that the X-ray emission originates close to the 
central black hole. In a number of Seyfert galaxies, the 2-20 keV 
continuum spectrum steepens with increasing X-ray flux (Iwasawa 
et al. 1996; Lamer et al. 2000; Lamer et al. 2003; Lee et al. 2000; 
Chiang et al. 2000; Done et al. 2000). This softening of the X-ray 
emission is commonly attributed to a change in shape of a single 
continuum component, possibly due to the Compton cooling of the 
X-ray emitting corona by the increased seed photon flux (Haardt 
et al., 1997). 

Recent spectral variability studies of MCG— 6-30-15 with 
RXTE (M c Hardy et al., 1998) and subsequently with ASCA (Shih 
et al., 2002) have shown that the photon index appears to saturate at 
high X-ray fluxes. To explain this effect, which Lamer et al. (2003) 
also observe in NGC 4051, M c Hardy et al. (1998) and Shih et al. 
(2002) independently proposed a two-component spectral model, 
consisting of a soft component with constant spectral slope but vari- 
able flux and a hard component with constant flux, so that at high 
fluxes the soft component dominates the spectrum and the spec- 
tral slope saturates to the slope of the soft component. This model 
naturally explains the observed spectral variability as a result of 
the relative changes in flux of the hard and soft continuum compo- 
nents, without recourse to any intrinsic change in the shape of either 
spectral component. Using time-resolved spectral analysis of a long 
XMM-Newton observation of MCG-6-30-15, Fabian & Vaughan 
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(2003) have shown that, on 10 ks time-scales, the X-ray spectral 
variability of this AGN is accounted for by a two-component model 
where the soft varying component is a power-law (with roughly 
constant slope) and the hard constant component is produced by 
very strong reflection from a relativistic disk. Fabian & Vaughan 
(2003) suggest that the large amplitude of reflection may be boosted 
by gravitational light-bending effects close to a Kerr black hole. 

As an alternative to the two-component model, Zdziarski et al. 
(2003) have suggested that the apparent saturation of the photon in- 
dex at high fluxes may instead be a result of a simple pivoting of 
the X-ray spectrum about some relatively high energy (~ 60 keV 
or greater), so that the photon index increases with the logarithm of 
the flux. In order to understand the physical implications of spec- 
tral variability in Seyfert galaxies, it is important to distinguish be- 
tween the spectral-pivoting and two-component models. Here, we 
present a model independent technique which enables us to de- 
termine conclusively that the two-component model for spectral 
variability is correct in the case of MCG— 6-30-15 and the Seyfert 
galaxy NGC 3516, but not in NGC 4051 which shows evidence 
for pivoting. This technique is outlined in Section 2. The data col- 
lection and reduction is described in Section 3. The observational 
results and the implications for the interpretations of spectral vari- 
ability in AGN are discussed in Sections 4 and 5. 



2 METHOD 

Consider the total fluxes measured in hard and soft energy bands 
(Fl, and F s respectively). We can express each flux as the sum of 
constant (C/, and C s ) and variable (//, and f s ) components. 
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Fg_ 5 (counts/s) 

Figure 1. Flux-Flux represen- 

tation of spectral variability as a result of spectral pivoting for a number of 
pivot energies (from the bottom E p = 30, 100,300, 1000,3000, 10000 keV 
with Fh vs F s relations described by power-laws with corresponding slopes 
a = 0.56,0.73,0.79,0.83,0.86,0.88). The flux-flux relations were simu- 
lated using the XSPEC FAKEIT command to create fake RXTE spectra, from 
which the fluxes in each band can be determined. 



Fs — fs "H C s 
Fh = fh + Ch 



(1) 
(2) 



A simple parameterisation of flux-dependent spectral variability re- 
lates the hard and soft variable fluxes, by //, = kff, where k and a 
are constants. If a < 1 this relation describes power-law pivoting 
spectral variability, where the varying power-law pivots about some 
energy F p (larger than the observed energies) while keeping the flux 
density at F p constant (see appendix in Zdziarski et al. 2003, and 
Fig. 1). As pivot energy increases, so a asymptotically approaches 
1. In the case of a = 1, the spectral shape of the varying compo- 
nent is constant with a hardness ratio k. Values of a > 1 correspond 
to a pivot energy which lies below the soft energy band. Incorpo- 
rating varying and constant continuum components, we obtain the 
following relation between F/, and F s : 



F h = k(F s -C s ) a + C h 



(3) 



By plotting vs F s (which we call a 'flux-flux plot'), we reveal the 
nature of the spectral variability. Specifically, if the varying com- 
ponent has constant spectral shape the Fh vs F s relationship will be 
linear with gradient k (i.e. the hardness ratio of the varying com- 
ponent), and an intercept C = Q, — kC s (this linear form occurs ir- 
respective of the actual spectral shape of the varying component). 
In contrast, if the spectral shape of the varying component depends 
on flux the Fh vs F s relationship will not be linear. In particular, 
if the spectral variability is explained by simple power-law piv- 
oting, with no constant components, the Fh vs F s relationship will 
be well fitted with a simple power-law. Other non-linear shapes of 
the Fh vs F s relationship might be envisaged for more complicated 
forms of spectral variability, or if the varying component spectrum 
is not a power-law. Note that such an approach has previously been 
used by Churazov et al. (2001) to show the existence of a constant 
soft spectral component in the high/soft state spectrum of the black 
hole X-ray binary Cyg X-l. 
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Figure 2. The flux-flux plot for MCG-6-30-15 (see text for details). The 
dashed line represents the linear fit to the data 



3 OBSERVATIONS AND DATA REDUCTION 

In order to test our method, we first use PCA data from RXTE ob- 
servations of MCG— 6-30-15 to examine the spectral variability of 
this source using flux-flux plots. MCG— 6-30-15 is an ideal test 
candidate as a two component model has already been suggested 
to explain the spectral variability of this source (M c Hardy et al. 
1998; Shih et al. 2002; Fabian & Vaughan 2003). We combine data 
from monitoring observations (Uttley et al., 2002) and a ~ 300 ks 
long-look observation (Lee et al., 1999), all obtained during RXTE 
gain epoch 3. We extracted 256s-binned light curves for the en- 
tire data set in two energy bands, 2-5 keV (soft, channels 0-13) and 
7-15 keV (hard, channels 19-41), using Standard 2 data from layer 
1 of PCUs 0,1,2 and standard Good Time Interval selection crite- 
ria (Lamer et al., 2003). The latest L7 background models for faint 
sources were used to determine the background level for the PCA. 



4 RESULTS 

4.1 Flux-Flux plots 

Using the 256s-binned hard and soft light curves, we made a flux- 
flux plot of the hard flux against the corresponding soft flux for 
each 256s segment (Fig. 2). A linear relationship between the soft 
and hard flux is apparent from visual examination of Figure 2. We 
next fitted a simple linear model, F/, = kF s + C, to the data. We ob- 
tained fit parameters of k = 0.538 and C = 1.170 but the fit is poor 
(acceptance probability <C 0.1%), due to significant scatter in the 
values of F/, for a given F s . However, the scatter is small and con- 
tains no systematic deviations from a linear model, implying some 
weak spectral variability which is uncorrelated with flux. 

The nature of the Fh vs F s relationship can be seen with 
greater clarity by binning up the Fh vs F s plot (Fig. 3). The lin- 
ear model then provides a good fit to the data (x 2 of 11.03 for 16 
degrees of freedom) for fit parameters of k = 0.5 16 ±0.006 and 
C= 1.312±0.038. We then fit the data with a simple power-law 
to determine whether the spectral variability may be caused by 
spectral pivoting of a single component in the abscence of any 
other spectral components. The fit is not formally acceptable (re- 
jected at P > 99.8% confidence). Fit parameters of both linear 
and power-law model fits are summarised in Table 1. If a con- 
stant spectral component is included in the power-law fit (i.e. we fit 
the data with Equation 3) the fit converges to a = 1, equivalent to 
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Table 1. Results of the linear and power-law fits to the flux-flux plots for each of the Seyfert galaxies in our sample. 
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Object 




Linear model 




Power-law model 








C 


X 2 /d.o.f. 


kpow 


a 


X 2 /d.o.f. 


MCG-6-30-15 


0.516±0.006 


1.3 10 ±0.036 


11.12/14 


1.251 ±0.030 


0.707 ±0.013 


34.94/14 


NGC3516 


0.598 ±0.008 


1.813±0.047 


13.45/14 


1.671 ±0.040 


0.662 ±0.014 


36.98/14 


NGC 5506 


0.822 ±0.022 


1.429±0.126 


10.62/14 


1.369 ±0.070 


0.850 ±0.020 


13.46/14 


NGC4051 


0.502 ±0.008 


0.452 ±0.022 


63.18/11 


1.040 ±0.003 


0.642 ±0.022 


22.95/11 



F 3 _ 5 (counts/ s) 

Figure 3. The binned flux-flux plot for MCG-6-30-15 . The dashed line 
represents the linear fit to the data and the dotted line the power-law fit 
(without a constant component). 



the two-component linear model. The 90% confidence lower limit 
for a using the power-law plus constant component model is 0.9. 
This value of a corresponds to a pivot energy E p > 10000 keV 
which seems physically implausible since the high-energy cut-off 
in MCG-6-30-15 is at ~ 160 keV (Guainazzi et al., 1999). We 
therefore reject the pivoting model for the spectral variability of 
MCG-6-30-15 . 

The linearity of the flux-flux relationship clearly implies a 
varying soft component with a spectral shape which is uncorrelated 
with flux (and which has a hardness ratio ~ 0.52). A positive inter- 
cept on the hard flux axis, suggests the existence of a constant-flux 
component, which is harder than the varying component. We note 
that the hard component may be weakly varying, in order to provide 
the observed scatter in the unbinned flux-flux plot or, alternatively, 
the soft component may show variations in shape which are uncor- 
related with flux. However, for simplicity we refer to these compo- 
nents as the soft varying component (SVC) and the hard constant 
component (HCC), which we identify with the components of the 
two-component models suggested by M c Hardy et al. (1998) and 
Shih et al. (2002), to explain the observed spectral variability of 
MCG— 6-30-15 . We will consider the specific spectral shapes of 
these components in Section 4.3. 

4.2 Flux-Flux plots for other Seyferts 

Having demonstrated that the two-component model for X-ray 
spectral variability is correct in the case of MCG— 6-30-15 , we 
now consider the spectral variability of some other Seyfert galax- 
ies. We show binned-up flux-flux plots for NGC 5506, NGC 3516 
and NGC 405 1 in Fig. 4, made using RXTE monitoring data, in the 
same manner as the MCG— 6-30-15 flux-flux plot (see Lamer et al. 



2000, Lamer et al. 2003 and Uttley et al. 2002 for descriptions of 
the data). We fitted all the flux-flux plots with the linear two compo- 
nent model and the simple power-law model. The fit parameters are 
shown in Table 1 and the corresponding models plotting in Fig 4. 
The NGC 35 16 flux-flux plot is well fitted by a linear model, but a 
simple power-law is rejected by the data (at P > 99.9% confidence). 
As in the case of MCG— 6-30-15 , the addition of a constant com- 
ponent to the power-law model yields a value of a = 1 (equivalent 
to the two-component model). At 90% confidence the lower limit 
for a is 0.86, which corresponds to a pivot energy E p > 3000 keV. 
Hence we rule out spectral pivoting in the case of NGC 3516. Like 
MCG-6-30-15 the spectral variability in NGC 3516 results from 
two spectral components. The NGC 5506 flux-flux plot can be sat- 
isfactorily fitted by both the linear and power-law models. This 
ambiguity probably stems from the narrower flux range covered 
by the NGC 5506 data (which is a factor of two lower than seen 
in MCG-6-30-15 and NGC 3516). We note that the best fitting 
power-law index to the NGC 5506 flux-flux plot corresponds to a 
pivot energy E p ~ 3000 keV which again we consider unlikely. 

A linear fit to the flux-flux plot for NGC 4051 can be re- 
jected. The power-law model is a much better fit but is still re- 
jected at P > 98%, although we note that the pivot energy predicted 
by this model (around 50-60 keV) coincides with that estimated 
from the photon-index-flux correlation by Zdziarski et al. 2002). 
The addition of a constant component to the power-law results in 
a formally acceptable fit (X 2 of 13.38 for 9 degrees of freedom), 
with fit parameters a = 0.738±g (corresponding to pivot en- 
ergy E p ~ 300 keV), k pow = 0.838 ± 0. 104 and C h = 0.183 ± 0.01 
for C s = 0. A hard constant component appears to be required but 
the bulk of the spectral variability of NGC 405 1 is a result of spec- 
tral pivoting. 



4.3 The nature of the hard constant component 

We have shown that there exist at least two modes of spectral vari- 
ability in Seyfert galaxies: variability described by a simple linear 
two-component model where the soft component (SVC) varies in 
amplitude but not spectral shape, with respect to a hard component 
(HCC); and spectral variability where the varying component piv- 
ots, but a hard component is also present. In either case the hard 
component may be constant or weakly varying, but its flux is not 
correlated with the flux of the soft component. The SVC may sim- 
ply understood as the varying power-law emission component (e.g. 
as measured using the 'difference spectrum', Fabian & Vaughan 
2003). An interesting question is the nature and physical origin of 
the HCC. We now briefly investigate the shape of the HCC in the 
case of MCG— 6-30-15 , using information contained in the flux- 
flux plot. 

As described in Section 2, the intercept on the hard flux axis 
of the flux-flux plot, C = Ch — kC s , where Q, and C s are the fluxes 
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Figure 4. Binned flux-flux plots for NGC 5506, NGC 3516 and NGC 4051. The dashed lines represent linear fits to the data and the dotted line the power-law 
fits (without a constant component except in the case of NGC 405 1 where we show the best fitting power-law plus constant model). 



of the HCC in hard and soft bands respectively, and k is the hard- 
ness ratio of the SVC. In principle, one could measure many values 
of C from flux-flux plots for many energy bands, to create a set of 
simultanous equations for the values of C which can be solved to 
yield the HCC flux in each band. However, the resulting equations 
(which we do not show here, but are simple to derive) are highly 
sensitive to the measured values of k and C so that even small un- 
certainties in the measured flux-flux relation make this method un- 
feasible. Instead, one can make a simple estimate of the shape of 
the HCC by making flux-flux plots in many bands and assuming the 
value of C s in the softest band and then extrapolating the observed 
linear flux-flux relation to find the value of hard flux corresponding 
to that C s . 

To implement this method, we measured lightcurves in 10 nar- 
row energy bands which do not overlap. The softest band covers the 
2-4 keV range, seven bands of width ~ 1 keV cover energies from 
4-11 keV and two further bands cover 11-13 keV and 13-15 keV. 
We then made flux-flux plots for each of the bands above 4 keV 
versus the 2-4 keV band. All the flux-flux plots were well fitted by 
linear relations with intercepts on the hard flux axis. We chose 6 test 
values of the 2-4 keV HCC flux, C s , in intervals of 0.3 count s _1 
from C s = count s _1 to C s = 1.5 count s _1 and determined the 
corresponding HCC flux in the other bands using the flux-flux rela- 
tions (errors are calculated by propagating through the uncertainties 
in the linear fit parameters). In Fig. 5 we show the resulting family 
of HCC spectra, normalised relative to the total average spectrum, 
in order to show the shape of the HCC independent of the instru- 
mental response. Note that the shape of the HCC above 4 keV is 
only weakly dependent on the assumed values of C s , so the method 
appears to be quite robust in constraining the HCC shape. It can be 
seen that the HCC is significantly harder than the total spectrum. 
Furthermore, a prominent feature can be seen around 6-7 keV, sug- 
gesting that the HCC contains stronger iron features than are seen 
in the total spectrum. 

A simple interpretation of the HCC is that it corresponds to 
reflection of the SVC from distant cold material (perhaps the puta- 
tive molecular torus), and hence can remain constant despite large 
changes in the illuminating continuum flux. However, we note that 
even the minimum possible contribution of the HCC (correspond- 
ing to C s = 0) exceeds 30% of the total flux at 10 keV, significantly 
larger than the < 10% of total flux expected for reflection from 
cold material which covers half the sky as seen from the contin- 
uum source (see Fig. 5 and George & Fabian 1991; Magdziarz & 



C B = 1 -5 
C a =1.2 
C 5 =0.9 
C s =0.6 
C s =0.3 
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Figure 5. Estimated HCC spectra normalised by the total average spectrum 
(see 4.3 for details). Also plotted as a stepped line is the expected maximal 
reflection spectrum (including a narrow iron Ka line) for cold distant mate- 
rial which covers half the sky as seen from the continuum source (assuming 
solar abundances and reflector inclination angle of 0°). 



Zdziarski 1995). We will examine the shape of the HCC in more 
detail, and in other AGN, in a future paper, but for now we con- 
cur with the interpretation of Fabian & Vaughan (2003) that the 
HCC may represent reflection from the inner disk which is boosted 
by light-bending effects close to a rapidly rotating black hole (al- 
though how such a component remains constant is not yet clear, but 
see Fabian & Vaughan 2003 for further discussion). 



5 DISCUSSION 

We have introduced the flux-flux plot as a model independent tech- 
nique for distiguishing between two modes of spectral variability in 
Seyfert galaxy spectra; pivoting of a single spectral component or 
changes in the relative normalisation of a constant hard and varying 
soft component with constant spectral shape (the two-component 
model). We have shown that the two-component model is correct 
for MCG-6-30-15 and NGC 3516 but the pivoting model applies 
in NGC 4051 (albeit with an additional constant hard component). 
The data for NGC 5506 is less clear cut and does not allow a defini- 
tive conclusion as to the nature of the spectral variability in this 
AGN, probably due to the smaller flux range sampled by the data. 
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This ambiguity in the nature of the flux-flux plot can be avoided by 
comparing fluxes in more widely seperated hard and soft bands, so 
that the index of the power-law flux-flux plot expected from pivot- 
ing is decreased (due to the greater relative change in the hard and 
soft fluxes). XMM-Newton observations should be able to provide 
the necessary separation of the hard and soft energy bands that is 
required, to distinguish the two possibilities in NGC 5506 and sim- 
ilar cases. 

It is interesting that both pivoting and constant-shape modes 
of spectral variability are seen in the black hole candidate Cyg X- 1 
(Zdziarski et al., 2002) with pivoting being observed in the hard 
state only and a constant spectral shape seen in both hard and soft 
states. Zdziarski et al. (2002) review how these different types of 
spectral variability can result from a two-phase disk-corona model, 
where the accretion disk illuminates the hot corona with low energy 
seed photons to produce the hard X-ray emission via inverse Comp- 
ton scattering (e.g. Haardt & Maraschi 1993), which in turn feeds 
back to affect the seed photon emission by heating the disk. For ex- 
ample, spectral pivoting results when the luminosity of the corona 
(which dissipates most of the gravitational energy release) is kept 
constant, while the seed photon luminosity varies. Spectral pivoting 
may also be produced if the corona is pair dominated. A constant 
spectral shape (as in the case of the two-component model) is pro- 
duced when the total luminosity of both disk and corona increases. 
Zdziarski et al. (2002) suggest that changes in the seed photon lu- 
minosity or total luminosity may correspond, respectively, to a re- 
duction in disk inner radius (so the corona sees more of the disk 
emission) or an increase in accretion rate. However, it is not clear 
how these physical models can explain the very rapid spectral varia- 
tions in AGN, where most of the flux variations (and corresponding 
spectral changes) occur on short time-scales (hours to days), much 
shorter than the expected viscous time-scales in these objects. 

Applying the AGN/X-ray binary analogy suggested by the re- 
markably similar variability properties of these objects (e.g. Ut- 
tley et al. 2002), one might think that the spectral pivoting ob- 
served for NGC 4051 implies that this AGN occupies the hard 
state, corresponding to a low accretion rate (few per cent of Ed- 
dington). However, the X-ray timing properties (power spectrum) 
of NGC 405 1 show a strong similarity to the soft state of Cyg X- 
1 (M c Hardy et al., 2003, in prep), and the low mass of this AGN 
suggests a high accretion rate of a few tens of per cent or larger 
(Shemmer et al. 2003, M c Hardy et al., 2003, in prep.). The power 
spectrum and spectral-timing properties (coherence, phase lags) of 
NGC 405 1 M c Hardy et al. (2003, in prep.), are also similar to those 
of MCG-6-30-15 (Vaughan et al., 2003), despite the apprently dif- 
ferent spectral variability modes operating in these AGN. 

One ap- 
parent difference between NGC 4051 and MCG-6-30-15 is that 
NGC 4051 shows a significantly larger variability amplitude (by 
about a factor 2) than MCG— 6-30-15 , on all time-scales. It is dif- 
ficult to see how such a large variability amplitude (especially on 
short time-scales) could be caused by variations in the seed pho- 
ton flux. However, if the seed photon luminosity is dominated by 
internal viscous heating in the disk (i.e. X-ray heating is not signif- 
icant) but the corona dissipates most of the accretion power, then 
a relatively small change in dissipation rate in the corona, if fed 
back to the disk, could result in a large change in the seed photon 
luminosity and relatively little change in coronal (and total) lumi- 
nosity, thus satisfying the conditions for spectral pivoting. Different 
dissipation mechanisms, and different coupling strengths between 
disk and corona could then result in the different spectral variability 
modes that we see. The similar power spectral shape and spectral- 



timing properties observed in MCG— 6-30-15 and NGC 4051 may 
indicate that these properties are linked to the more general proper- 
ties of the corona (e.g. radial distribution of emission, temperature 
gradient) rather than the disk-corona coupling mechanism (e.g. see 
Kotov et al. 2001 for a possible model which can explain these tim- 
ing properties). 

Finally, we note that the hard, constant component in the spec- 
trum of MCG— 6-30-15 appears to be broadly consistent with a 
pure reflection spectrum but the strength of this component sug- 
gests that the reflection is somehow boosted, perhaps by relativistic 
beaming and gravitational light bending effects (Martocchia et al., 
2000). Such a scenario can be realised if the reflector lies close to 
a central Kerr black hole, as also suggested by Fabian & Vaughan 
(2003) who describe how such a model might explain the constancy 
of the reflection spectrum despite large changes in the soft compo- 
nent flux. 
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